Low Temperature Transport Properties of Strongly Interacting Systems 
- Thermal Conductivity of Spin-1/2 Chains 
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We outline a general approach to the computation of transport properties of interacting systems 
at low temperetures and frequencies. We show that if the fixed point and the irrelevant opera- 
tors around it are known, then by studying the structure of the softly violated conserved currents 
chracterizing the fixed point one may set up an effective calculation in terms of a memory matrix 
formalism. We apply this approach to the computation of thermal conductivity of spin chains em- 
bedded in a matter matrix and interacting with its phonons. The results are found to be in very 
good agreement with experiment. 



The study of transport properties of strongly interacting 
systems has been of great theortical and experimental 
interest for a long time. We shall concentrate in this 
contribution on developing an effective approach to the 
problem when the system is not far from its fixed point; 
in other words, when it is probed at low temperatures 
and at low frequencies. Subsequently we shall carry it 
out in detail for a system of spin 1/2 chains coupled to 
phonons. 

Our approach consists of the following elements: 

1. Identify the fixed point of the hamiltonian which de- 
scribes the system, as well as the irrelevant operators 
around it. 

The fixed point itself is typically insufficient to describe 
low energy transport properties; it is scale invariant and 
translationally invariant and thus unable to degrade a 
current, leading to infinite conductivity. To obtain finite 
conductivity one needs to take into account terms which 
break translational invariance and, more generally, vio- 
late the conservation laws associated with the fixed point. 

2. Study the (weakly violated) conserved charges around 
the fixed point hamiltonian. 

A fixed point H* is scale invariant and often has several 
conserved quantities P, [P, if*] = 0, associated with it. 
In 1-d, for example, if the fixed point is conformally in- 
variant it has an infinite number of conserved quantities. 
When the irrelevant operators around the fixed point are 
taken into account, most of these quantities no longer 
commute with the low energy hamiltonian. The conser- 
vation of most of these quantites is strongly violated, but 
some may be only weakly violated and then significantly 
influence the low energy dynamics of the system. 
Typically the current whose correlations determining the 
transport properties under considerations will be among 
those almost conserved charges, or "protected" by them 
in the following sense. When a system possesses some 
conserved quantities P, these may "protect" the cur- 



rent J from degrading (this occurs when the cross- 
susceptibility xjp 7^ 0) leading to a pure (i.e. S(u)) 
Drude peak and infinite d.c. conductivity. When the 
conservation of the pseudo-momenta P is softly violated 
they will, instead, lead to very long time tails in the de- 
cay of the current J. This occurs since states with a 
finite pseudo-momentum P typically carry also a finite 
current J since \JP 0- The component of the current 
"parallel" to P, J\\ P = (xpj/xpp)P W1 U therefore decay 
slowly. The presence of such approximately conserved 
quantities leads then to a natural hydrodynamic descrip- 
tion of the system where a separation of fast and slowly 
decaying modes takes place and a consistent scheme of 
calculation of the slow mode conductivities can be carried 
out in terms of matrix of decay rates of these modes. 
Let us explain and illustrate these ideas in more detail. 
We begin by arguing that given a conserved charge P, it 
will "protect" J if xpj 7^ 0. Indeed, imagine preparing 
at t = a state carrying a current (J). Then necessarily 
that state will also have a non vanishing (P), 

{P) = —(J). 
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FIG. 1 

In the limit t — > oo since P is conserved its expectation 
value will not have changed, and it will in its turn induce 



a non vanishing expectation value for </, 

-(J(t = 0)) 



Y 2 
AJP 



lim(J) = —{P) = 
t^oo Xpp XppXJJ 



(see Fig. 1). Since the current tends asymptotically to 
a constant value we find that the conductivity will have 



a (5(a>)-Drude peak containing a fraction 



XppXJJ 



of the 



1 y 

total weight; in other words, a Drude weight D = ^^j- 
An immediate consequence is that integrable models, 
having infinitely many conserved quantities, will typically 
also have an infinite dc conductivity. 
If the charge P is not conserved but slowly decaying it 
will induce slow (long-time) decay in J (see Fig. 2). 
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where n s is the spin contribution obtained after subtract- 
ing the phonon background from k^. The relevant energy 
scales of the system are: The Debye temperature char- 
acterizing the phonons - 8^ ~ 400 K and the spinon 
interaction scale J/ kp ~ 2600 K characterizing the spin 
chain. 




FIG. 3 



FIG. 2 

This slow decay will show up as a peak in (j{lo) and the 
decay-rate Tp of P will determine the width of the peak 
and the dc value <r(0) « D/Tp of the optical conductiv- 

2 

ity, the area under the peak being again, P> = \ • 
3. Having identified the slowly decaying charges of the ef- 
fective model describing the low energy properties of the 
system we can compute the low temperature dc trans- 
port using a method that separates fast and slow modes, 
incorporating the former in the dynamics of the latter. 
Such a method is the memory matrix approach which can 
be used very efficiently and controllably when combined 
with the RG considerations outlined above. 

Let us apply this approach to a system of spin chains 
embedded in some 3-d lattice and interacting with its 
phonons. Such systems (including, in particular, various 
compounds of SrCuO) have been recently studied in de- 
tail by, e.g., Sologubcnko et al^ The authors have mea- 
sured the heat conductivity along the three main axes 
of the sample. They observed that while conductivities 
along the a and c axes almost coincide, the conductivity 
along 6, the axis along which the spin chains lie, has an 
enhancement which they interpret as being due to contri- 
bution of the spin degrees of freedom (Fig. 3, for exam- 
ple, presents data from a corresponding measurement in 
Sr 2 Cu0 3 ). In the temperature range 60 K < T < 200 K, 
they gave the fit: 

k s (T) ~ exp(TVT) , 
T* w 0.429Z5 



One immediate question that arises is why is k s actu- 
ally determined by the lower scale Qr>- Furthermore, by 
plotting the above data of In k s and of the pure phonon 
contribution lnK a c vs. 1/T, one observes that the slope 
of the latter is larger by a factor of 2. Again, why? 
We shall find that a rather subtle interplay of (approxi- 
mate) conservation laws and quantum dynamics under- 
lies the experimentally observed heat conductivity, and 
the approach outlined above is necessary to fully account 
for it. 

We begin by discussing the low energy effective hamilto- 
nian. First consider a single spin chain, 

N N 

H * = 2 E ^ ( s i~ s i + s ^ s t) + E j 'j s * s ! 

As is well known, 2 for spin chains with short range in- 
teractions the fixed point hamiltonian is the Luttinger 
liquid, 

H L l = -i(Ja) J dx(^ R d x tf} R - ipld x tf} L ) + J z J dxp(x) 2 

with tpR./L being right/left moving fermi field, and J and 
J z are some average values of and J?- respectively. As 
a reminder, the fixed point can be obtained by carrying 
out a Wigner- Jordan tranformation and then linearizing 
the resulting fermions around the Fermi points ±fcp (see 
Fig. 4). 















left mover 


right mover 



with q the lattice deformation parallel to the spin chains 
direction, and P the conjugate momentum. 
Integrating the corresponding propagator over the per- 
pendicular directions, we obtain the propagator along the 
chains: 



FIG. 4 

It will be convenient to change to bosonic variables in 
terms of which the fixed point hamiltonian takes the form 
(to leading order in \J Z \/J), 



H 



LL 



where 



J 



K 



1 



2J, 
nJ 



We now consider the irrelevant operators around the 
fixed point. The irrelevant operators come from vari- 
ous sources. We shall divide them to Umklapp and non- 
umklapp operators. To the latter category belong local 
terms coming from band curvature around kp: for exam- 
ple J ip' R d 2 ipii. However, the operators that are impor- 
tant for transport are the Umklapp operators that reflect 
the underlying lattice structure. Only they break trans- 
lation invariance and can degrade the currents. They 
have the following structure: 



u 

nrn 



9nm / dx[e iAk ™ x Jl il>U x + 3o-)^l{x + ja) + h, 

3=0 

- ( dx[e iAknmX e i2n ^ xS > + h.c] 



~~ (2na) 

where A.k nm — n2kp — mG is the momentum transfer as- 
sociated with the Umklapp process where n particles are 
transferred from one Fermi point to another while giving 
up to the lattice m units of lattice momentum G = 2ir/a. 
The particular values of the couplings g% m depend on the 
couplings in the microscopic model. These terms are ir- 
relevant perturbations around the fixed point not only 
by power counting but also because the x-dependent ex- 
ponentials in them suppress their contribution exponen- 
tially. However they are the only source of dissipation. 
We now consider the complete system consisting of an ar- 
ray of parallel spin chains interacting with 3-dimcnsional 
acoustic phonons. The 3-d phonon system projected 
along the axis describing deformations of the lattice par- 
allel to the chains is described by 



Hl D = 



d 3 x 
~27 



(ttP) 2 



d 2 k. 



(« 2 + E M «2^) 



ln[(u; 2 



In real space it takes the form l/(x 2 + v 2 t 2 ) with v p 
being the slowest phonon velocity. This again is the 
propagator of a K = 1 Luttinger liquid describing the 
phonons, H p — v p J || ((tHI) 2 + (d x q) 2 ) . Again we 
need to add to the combined phonon-spinon fixed point 
H* = J2 a + ^2 a H" L (the summation a is over the 
spin chains) all the irrelevant operators. As before they 
fall into two categories, the Umklapp and non Umklapp 
operators. The most important ones of the former one 
have the form, 



H 



U,s—p 



q u, P 



dx[t 



iAk n 



h.c 



(2na) n 

(where <fi is the bosonic field in a particular chain), while 

ion" 



an example of the latter is H™° nU = J (dcb) 2 d x q. 



We now turn to our main interest, computing the ther- 
mal conductivity k s (lj,T) of the spin chains coupled to 
phonons. The thermal conductivity can be expressed in 
terms of a heat current correlation function, 

K (u,T) = (J Q , Jq)(u,T)/uj 

where Jq, the heat current expressed in bosonic vari- 
ables, is J Q = -J2al dx v 2 U a d x cb a - J d 3 xv%Pd x q . 
The correlation is to be computed with respect to the 
I low-energy hamiltonian, 



H, 



low — E 



= H* 



H 



nonU 



As explained in the outline we need to identify the (ap- 
proximately) conserved "charges" of the low-E Hamilto- 
nian, to find out whether they induce a slow decay of the 
heat current. We now show that the quantities, 

J s = vK^T dx [ip f Ra ^R a - ^ La 1pLa] = / dxll c 

a J ct 

P T = - £ / dxll a d x (f> a - / d 3 xPd x q 

where, J s is the spin current and Pt the momentum op- 
erator, are the "slow modes" which in turn protect Jq 
rendering it slow too. 

Indeed, J s and Pt commute with Hll and with H nonU 
and their conservation is violated only through H u with 
which they do not commute, thus inducing a slow current 
decay. More importantly, certain linear combinations of 
J s and Pr, the "pseudo-momenta" 



Pn 



1 

2^i 



^ kn . m Js + Pt 



decay even slower as they commute with Hll + H nonU + 
H^ m + H^£-p and are therefore exactly conserved if 
only a single type of Umklapp with quantum numbers 
n and m is present. We note that the pseudo momenta 
can be written as P nm — P\ a t + ^G{Nr — Nl) with 

Plat = Efc fc 4 c fc « k F (N R - N L )+P T being the lat- 
tice momentum. Unlike other (approximately) conserved 
quantities, they decay exponentially slowly with the tem- 
perature as their violation requires processes away from 
the Fermi energy £ 

The heat current on the other hand does not commute 
with both H nonU and H u and would therefore decay fast, 
but is protected by J s , Pt and their linear combinations 
since XJ q ,J s ,XJ q ,p t 0. 

We proceed to discuss transport in the presence of several 
approximately conserved - "slow" - variables: Ji, J2---Jn- 
We shall introduce the memory matrix formalism^ which 
is very effective under the circumstances as it allows the 
separation of the slow modes (J s , Py, Jq, in our particu- 
lar case) from the fast modes. 

To set up the formalism one introduces a scalar product 
in the space of operators of the theory 

(A(t)\B) = ±J*d\(A{t?B(i\)). 

In terms of this scalar product one can express the dy- 
namic correlation functions as follows, 



A 



iT 



dte iwt (A(t)\B) 



B 



uj — C 



dte wt ([A(t),B]) 



(A\B) 



where the Liouville operator £ is defined as LA = [if, A] . 
The matrix of conductivities (Kubo formula) is then, 



a pq (u,T) = —C. Jp . Jq {uj) 
(p, q = 1 • • • N) . In our case the thermal conductivity is 



K S (U!,T) = -<Jqq{uJ,T). 

However the dc conductivity has no good perturbative 
expansion: a ~ with T the decay rate of the current 
is singular in perturbation theory. Put in other words, if 
we wish to compute the conductivity as a perturbative 
expansion of the irrelevant operators around the fixed 
point, this would be an arduous task in view of the fact 
that the conductivity computed from the fixed point is 
infinite. Furthermore, even the perturbative expansion 
for l/a(uj — 0) turns out to be ill behaved in the pres- 
ence of slow modes. While the short-time decay rate (see 
Fig. 2) of the current is perturbative, the oj — > limit 



requires to take into account the presence of approximate 
conservation laws. 

We seek therefore to compute a quantity which has a 
good perturbative expansion, and to this purpose intro- 
duce M(cv,T) - the Memory Matrix, essentially the ma- 
trix of relaxation rates. 
The matrix is defined as 



Q 



d t Jq 



where Q is the projection away from slow modes 

Q = i-Y J \j P )hr 1 ) Pq (j q \- 



In terms of the memory matrix the conductivity matrix 
is, 

&(u, T) = X (T) (m(«, T) - iu>x{T)) ^ X(T) 
with x the susceptibility matrix, 



Applying the formalism in our case we find that the mem- 
ory matrix is a sum over the Umklapp processes (nm), 
given - to leading order in c^ nn - by, 



M 



1 

T 



-p) 



where (the matrix indices p, q take the values s, T, Q) 



M pq 



{FP;F"Z-(FP;FiZ =0 



M pq 

nm,s—p 



(ipP . F q VO _ (JpP . F q \0 
\ s—pi s—p/uj \ s—pi s—p/uj—0 



here F p = i[J p , H u ], F s % = i[J p , H U ' S ~ P ] and (F p ; 
is the retarded correlation function calculated with re- 
spect to H* . As the perturbative expansion is in irrela- 
vant operators with respect to the fixed point the expan- 
sion is expected to be rapidly converging at low temper- 
atures if slow and fast time scales are well seperated and 
all of the slowest modes have been taken into account. 
Carrying out the computation of the various correlation 
functions (see Ref. for details) we have: 



k,(T) w v 2 T 3 (M- 1 ) TT + 2(M- 1 ) QT + (M- 1 ) QQ 



with the typical matrix elements, 



~ [&k nm r 2K - 2) e- vAk - /2T 



M pq 

rvm,s—p 



rp(2n 2 K-l) e -v p Ak nm /2T 



Note that the spinon processes decay exponentially 
fast, with the exponent —vAk nm /2T, while the spinon- 
phonon exponents contain the much slower phonon ve- 
locity v p v, leading to a much slower decay with the 
exponent —v p Ak nm /2T. The latter will therefore clearly 
determine the thermal conductivity. 

But which of the scattering processes (n, m) will domi- 
nate? 

At low-T exponential factor prevails hence the smallest 
Ak nm . At this point our discussion must distinguish 
between commensurate and incommensurate magnetiza- 
tion, or in fcrmionic language commensurate and incom- 
mensurate filling. The magnetization can be in principle 
tuned by varying an external magnetic field h. 
Close to commensurate filling kp ss G^? 2 -, and then 
the dominant processes appear to be H^ omg where 
Afc n0i „ lo w 0. However, because of the conservation 
laws discussed earlier, it is the next leading term H% imi 
with Ak„ 1<mi = ±G/tiq which determines the decay rate. 
Technically, this arises as one of the eigenvalues of the 
matrix M is not affected by H% omo (as [P„ o m , H% omo ] = 
0) but determined by H% im . This smallest eigenvalue 
will then determine the size of M _1 and therefore the 
heat conductivity. The particular case of half filling is 
discussed in detail below. 

On the other hand, at a typical incommensurate filling 
it will depend on the temperature which processes are 
dominant and subdominant and we need to sum over all 
terms (do saddle-point approximation with respect to ni) 
and find, 

Atypical ~ exp[ C (/%G) 2/3 ] . 

with c a constant of order 1. 

We now turn to the experiment by Sologubenko et al. 
discussed earlier. The experiment was carried out at 
h = corresponding to half filling, Uf = 7r/2a = G/4. 
Therefore A&21 = (n=2, m=l). But we need at 
least two Umplapp terms and the next smallest is Ak = 
G/2,(n — 1, m = 0). Recall also that as v p <C v, it 
follows that M nm;S _ p 3> M nm . Hence, the dominant 
contribution to the thermal conductivity comes from, 
{M-^tt ~ l/M^ l s _ p {G/2,T) and we have 

k(T) » /t (T/T*) a ( 1 -^e T "/ T 

with 

T* = v p G/A . 

We find therefore that the second strongest rate wins 
(cf. the expression for M): it is determined by v p via 
a phonon process and is characterized by the momen- 
tum G/2. The G/2 transfer momentum characterizes 
the dominant spinon - phonon Umklapp process, and 



clearly distiguishes it from pure phonon Umklapp pro- 
cesses characterized by momentum transfer G. We ex- 
pect therefore that in the pure phonon thermal conduc- 
tivity (axes a and c in Fig. 3) the scale 2T* would appear. 

To compare our findings with the experimental data we 
need to express our expression in terms of ®d- As- 
suming an isotropic phonon dispersion one has: 8d w 
u„(67r 2 /a 3 ) 1/3 « 0.6 v p G. 
Therefore: 

T* w 0.46/5, (theory) 

T* w 0.426D, (experiment). 

Taking into account possible ambiguities in the fits to the 
experiments and that the phonon dispersion is probably 
not completely isotropic, part of this excellent agreement 
may be accidential. But further confirmation of our the- 
ory comes from the observation that the ratio of slopes 
of the spinon contribution compared to the pure phonon 
contribution (on a semilogarithmic graph of n vs. 1/T - 
see Fig. 5& ) is approximately 1:2 as discussed earlier. 
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We may also consider the effect of a magnetic field on the 
thermal conductivity. (To observe it experimentally one 
needs a material with much smaller spin energy scales J 
than SrCuO). The magnetic field modifies the value of 
&f according to 

fc F = ^(l + M)«^(l + V(^))- 

As the field h is varied the system passes through com- 
mensurate fillings, Ak nrn = n2kp — mG = and in- 
commensurate fillings, Ak nm = n2kp — mG ^ 0. Thus 
different Umklapp operators become effective leading to 
a fractal- like dependence on M « h/irJ (see Fig. 6). 
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FIG. 6 



We expect some dips to be experimentally observable. 
In conclusion, 

• Transport is strongly influenced by conserved 
"charges" : low energy processes cannot relax heat 
current. 

• Exponents are determined by slowest mode in the 
system: typically phonons. 

• Memory Matrix approach, separating slow and fast 
modes, allows controllable calculations. 

• Calculations fit experiments 

• Interesting predictions on the magnetization de- 
pendence of the heat transport have been made. 
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